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ENERGY TRANSFER AND TRIADIC INTERACTIONS IN COMPRESSIBLE 

TURBULENCE* * * § 


F. BATAILLE*, YE ZHOU*, AND .JEAN-PIERRE BERTOGLIO§ 

Abstract. Using a two-point closure theory, the Eddy-Damped-Quasi-Normal-Markovian (EDQNM) 
approximation, we have investigated the energy transfer process and triadic interactions of compressible 
turbulence. In order to analyze the compressible inode directly, the Helmholtz decomposition is used. The 
following issues were addressed: (I) What is the mechanism of energy exchange between the solenoidal and 
compressible modes, and (2) Is there an energy cascade in the compressible energy transfer process? It is 
concluded that the compressible energy is transferred locally from the solenoidal part to the compressible 
part. It is also found that there is an energy cascade of the compressible mode for high turbulent Mach 
number (M t > 0.5). Since we assume that the compressibility is weak, the magnitude of the compressible 
(radiative or cascade) transfer is much smaller than that of solenoidal cascade. These results are further 
confirmed by studying the triadic energy transfer function, the most fundamental building block of the energy 
transfer. 

Key words. Compressible turbulence, energy transfer, triadic interactions, closure theories. 

Subject classification. Fluid Mechanics 

1. Introduction. It is well known that compressible turbulence plays a prominent role in a wide range 
of important scientific and engineering applications, including high speed transport, supersonic combustions, 
and acoustics. Recently, a large body of publications have been devoted to study various aspects of com- 
pressible turbulence using both direct numerical simulations (DNS) and large eddy simulations (LES). Here 
we simply mention a few representative works, such as Feiereisen et al. (1981), Passot and Pouquet (1985), 
Lee et al. (1991), Sarkar et al. (1991), Erlebacher et al. (1992), Kida and Orszag (1990), Blaisdell et al. 
(1993) and Porter et al. (1994). For a comprehensive review, the reader is referred to Lele (1995). These 
numerical simulations have substantially improved our understanding of compressible turbulence. Never- 
theless, some basic physical processes of compressible turbulence, such as the energy transfer and triadic 
interactions, have not been explored even at low turbulent Mach number. For example, do we expect an 
energy cascade process of compressible velocity modes? How does the energy exchange process between the 
solenoidal and compressible modes take place? These type of studies require a substantial spectral scale 
range of interactions. As a result, it is very hard to utilize the DNS databases since these simulations are 
limited to very low Reynolds numbers and have only very limited spectral ranges. While LES can provide 

* This research was supported by the National Aeronautics and Space Administration under NASA Contract No. NAS1- 
19480 while the first and second authors were in residence at the Institute for Computer Applications in Science and Engineering 
(ICASE), NASA Langley Research Center, Hampton, VA 23681-0001. The first author was also partially supported by the 
Foreign Ministry of the French government. 

UNSA, Centre de Thermique de Lyon, UPRES A CNRS 5008, 20 av. A. Einstein, 69620 Villeurbanne, France (email: 
bataille @cethil. insa-lyon.fr ) . 

* Institute for Computer Applications in Science and Engineering, NASA Langley Research Center, Hampton, VA 23681 
and IBM Research Division, T.J. Watson Research Center, P.O. Box, 218, Yorktown Heights, NY 10598 (email: zhou@tcase.edu 
and/or yzhou@xuatson.ibm.com). 

§ Laboratoire de Mecanique des Fluides et d’Acoustique, URA CNRS 263, 36, av. Guy de Collongue, 69130 ECULLY, 
FRANCE 


1 



databases at substantially higher Reynolds numbers, subgrid models must be introduced. 

Another way to generate high Reynolds number databases is by using two-point closure models. Di- 
rect Interaction Approximation (DIA) of Kraichnan (1959) is a well-established approach. Many authors 
have studied incompressible turbulence using DIA (see Leslie (1973)). The method of Eddy-Damped-Quasi- 
Normal-Markovian (EDQNM) models (Orszag (1970)) has been shown as a simpler, but effective, alternative 
to DIA. Recently, Bertoglio et, al. (1996) have presented DIA and EDQNM equations for a weakly compress- 
ible turbulence. 

The aim of this paper is to use the EDQNM closure theory to study the energy transfer and triadic 
interactions of compressible turbulence. The paper is organized as follows: First, we review closure assump- 
tions of EDQNM and present the resulting transport equations. Second, we perform a detailed analysis of 
the non-linear transfer terms. Finally, we investigate the most fundamental aspect of the energy transfer 
process, the incompressible and compressible triadic interactions. 


• E ss 
A E tx (M = 1 ) 
l ' E (X '<M=10 l ) 


E <X (M ( =10' 2 ) 



HI 1 10* 10 1 HI 4 

Figure l:Solenoidal and compressible velocity correlations 
for different turbulent Mach numbers 


Fig. 1. 

2. Compressible EDQNM Model. The basic set of equations are the Navier-Stokes and continuity 
equations. The fluid is assumed to be homogeneous, isotropic and barotropic. Reynolds average and Fourier 
transform are used to obtain the fluctuating turbulent field equations in spectral space. The equations are 
partially linearized with respect to the density fluctuation which leads to the condition M t < 1, where M t is 
the turbulent Mach number defined a sM t — q 2 is twice the turbulent kinetic energy, and c 0 is the 

sound speed. To analyze the compressibility effects, we use the Helmholtz decomposition to split the velocity 
vector into a solenoidal part u s (K,t) } which corresponds to the velocity fluctuations perpendicular to the 
wave vector K in the Fourier space, and a compressible part u c (K,t), which corresponds to fluctuations in 
the direction of the wave vector. 

The classical DIA approach (two point - two times) is used to derive compressible DIA equations. These 
equations are the starting point used to derive the compressible EDQNM equations (in two points and one 
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time). The methodology is the same as that of the incompressible case (Leslie (1973) and Lesieur (1987)) 
but more equations are now needed. For more details, see Bertoglio et al. (1996) or Bataille (1994). 

2.1. Spectral Equations. The EDQNM governing equations of weakly compressible turbulence are 
the following: 

- an equation for the spectrum ( E ss ), the auto-correlation of the solenoidal part of the velocity field: 

(1) ^ = ~ 2vK2eSS ^ K ' *) + T ss {K, f); 

- an equation for the spectrum ( E cc ), the auto-correlation of the “purely compressible” part of the 
velocity field: 

(2) £ E CC (K, t) = -2 u'K 2 E cc (K, t) + T CC (K, t ) - 2c 0 KE CP (K , t); 

- an equation for the spectrum of the potential energy ( E pp ) associated with the pressure: 

(3) £ E PP (K, t) = 2c 0 KE CP {K, t); 

- an equation for the spectrum of the pressure- velocity correlation ( E CP ): 


(4) E cp (K, t) = -v'K 2 E CP (K,t)+T CP (K,t) + c 0 K(E cc (K } t) - E pp (K,t)). 

at 


In the case of a Stokesian fluid: 


( 5 ) 


, A + 2 fi 4 

V = = 77 V 

< p > 3 


fi and A, two dynamic viscosities, are assumed to be uniform. 

2.2. Energy Transfer Terms. In equations (l)-(4), T ss ,T CC , and T CP are the transfer terms. They 
contain several contributions: 


( 6 ) T ss = Tf s + T$ s + Tf s + Ti s + T* s 

(7) T cc = if c + T? c + T 3 cc + T 4 cc + Tf c + Tf c 

(8) T CP = Tf p + Tf p + Tf p + 7f p + Tf p + 7f p 

Different contributions appearing in the transfer term acting on the solenoidal field are: 

(9) Tf s = f e S K S p -ss-ss E s S{p t)E ss {Q t)dpdQ 

( 10 ) T 2 SS = [ £ 0 ff p -ss-cc E ss {p>t)E cc m)dpdQ 

(11) T 3 ss = - / ^ (x;y + z 3 )O s k % ss - ss E ss (K, t)E ss (Q , ()<£PdQ 

JA V 
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( 12 ) 


rf s = (2 xy)eff p - ( * s - cc E ss (K,t)E cc (Q,t)dPdQ 

(13) t£ s = - J ^ (z(i - 2 2 ))^pQ C ' c “- ss £ ss (^,t)£ ss (<3^)^Q 

Different contributions to the transfer term in the E cc equation are: 

(14) Tf c (K, t) = J^^({x + yz) 2 )6 c K c p - Q ss ~ ss E ss (P , t)E ss (Q, t)dPdQ 

(15) T.2 C (K, t) = J ^ 9 C K % SS - CC E SS {P, t)E cc (Q, t)dPdQ 

(16) T.£ c (K, = (x 2 )0%% cc - cc E cc (P, t)E cc (Q, t)dPdQ 

(17) Tf c (K, t) = - f ^ 2z(l - z 2 )e c K % ss ~ ss E cc {K, t)E ss (Q, t)dPdQ 

J A Q 

(18) Tf c (K, t) = - f (2z 3 - z + xy)e° c P n CC - ss E ss (Q, t)E cc (K, t)dPdQ 

J A W 

(19) Tf c (K, t) = j ^ (2 xy)6% c P Q CC - cc E cc {K, t)E cc (Q, t)dPdQ 

And lastly, the contributions to the transfer term in the pressure velocity correlation are: 

(20) Tf p (K, t) = jJ?L (x+ 2 ^ )2 0 kpq SS ' SS £ S 5 (^ *)^ SS (0> OdPdQ 

(21) T 2 cp (tf , *) = J 2(l~- V ^) ' 0 k° Q SS ~ CC eSS ( P ’ 

(22) Tf p (K, t) = 0 k C P q CC ~ CC E cc {P, t)E cc (Q, t)dPdQ 

(23) Tf p (K, t) = - [ ^ z( 1 - z 2 )0 p £q SS - ss £ cc (K, f)£ 5S (Q, t)dPdQ 

Ja Q 

(24) T$ P {K, t) = - J ^ 2g3 ~ 2 Z+Xy 0 p K % cc - ss E ss (Q, t)E cc (K, t)dPdQ 

(25) T 6 CP (K, 0 = f ^ ( xy)e pc P Q CC ~ cc E cc (K , t)E cc (Q , *)dPd(? 
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Figure 2:SolenoidaJ transfer term 

Fig. 2. 

The integration in the P, Q plane extends over a domain such that K, P and Q form a triangle. The 
expressions x,y,z are standard coefficients associated with the geometry of the triad and they the cosines 
of the angles respectively opposite to K, P, Q. Furthermore, temporal functions (defined by 0 y s) are the 
decorrelation times deduced by integrating the DIA equations over time. Their expressions can be found in 
Bertoglio et al. (1996). 

The transfer terms appearing at the transport equations for the solenoidal and compressible inodes 
will be analyzed in detail in the next section. Since there is no transfer term in the equation for the 
pressure-pressure correlation, this equation characterizes the exchange between the solenoidal mode and the 
compressible mode (E cc + E PP ). Finally, the effect of T CP is to interchange energy between E cc and E pp 
(Bataille and Bertoglio (1993)). 

2.3. Energy Spectra. A force is applied to the large scales of the solenoidal velocity. Our analysis is 
carried out when both solenoidal and compressible modes have reached their asymptotic stationary states. 

In Figure 1, we present the spectra of both solenoidal and compressible components of the velocity 
correlation corresponding after the asymptotic state is reached. The Taylor micro-scale Reynolds number, 
R Vi is approximately 140. Here R e = | and e is the dissipation. In the figures, the dimensional unit of 

the correlation spectra is given as m 3 s~ 2 and the dimensional unit of wavenumber K is defined asm -1 . The 
energy associated with the purely compressible mode is found to vary as the square of the turbulent Mach 
number. At low and moderate M t , the spectrum of the compressible component shows a K~ n ^ 3 behavior 
in the inertial range (Bataille and Bertoglio (1993) and Bertoglio et al. (1996)). This behavior has been 
confirmed with Large Eddy Simulation by Bataille et al. (1996). 

In order to examine the dependence of the turbulent Mach number, it was allowed to vary from 10 -4 to 
1. Note, however, strictly speaking, our model is valid only for small M*. Nevertheless, computations were 
carried out up to M t = 1 in order to study the limiting behavior of the model. 
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Figure 3: Contributions of the soienoidal transfer term 


Flri. 3. 



Fig. 4. 


3. Study of the Transfer Terms. 

3.1. Soienoidal Transfer Term. In Figure 2, we observe that T ss has the usual shape observed in 
incompressible turbulence studies. Specifically, is negative in the large scales and positive for small 
scales. Physically this corresponds to the energy transfer from the large scales to the smaller ones. 

There are 17 contributors to the transfer function of compressible turbulence. This should be be com- 
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Fig. 5. 



Fig. 6. 


pared with its counterpart of incompressible turbulence: 

(26) T S SS (K, t) = T? S (K, t) + Ti s (K, t). 

Indeed, the two key contributors to T ss are Tf s and T-f s , which are the same terms that one finds in 
incompressible turbulence. These terms are usually called “input” and “output” terms. For our weakly 
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Fig. 7. 



Fig. 8. 


compressible turbulence, they are much more important contributors than the others (see Figure 3). Indeed, 
other contributors are negligible in comparison with Tf s and T§ s . We stress that the summation of these two 
terms is the net energy transfer of incompressible turbulence (solid line). To clearly illustrate the behaviors 
of the new terms, we plotted these “compressible” contributors (T<f 5 , Tf s and T^ s ) in an enlarged scale 
(Figure 4) in order to observe their relative magnitudes. We found that T 5 S5 is much larger than T^ s and 
Tf s , and consequently it makes the largest contribution to the compressible effects in T ss . An important 
feature of this term is that it is negative for all spectral space, indicating that the energy transfer has been 
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Fig. 9. 



Fig. 10. 

transferred from the solenoidal mode to the compressible mode. 

We observe in Figure 5 that the compressibility has negligible effect on the solenoidal transfer term 
T ss . The reason is that the dominant terms in T ss , e.g., the incompressible contributions (T? s and T 3 55 ), 
are independent of M t (Figure 6). On the other hand, the “compressible” contributions (which are smaller) 
depend strongly on the values of the turbulent Mach number. Figure 7 illustrates the dependences of T|f ‘ 5 and 


9 






Fig. 11. 


T$ s on the turbulent Mach number. We found that all the spectra of T$ s can be collapsed by dividing M 2 
(Figure 8). Therefore, we conclude that T§ s is proportional to M 2 , a result that can be found analytically. 
Since the dominant contributions in T ss are insensitive to the variations of turbulent Mach number, T ss is 
not affected by the compressibility. 

3.2. Compressible Transfer Term. In this subsection, we will study the transfer term T cc , and 
its individual contributors. This compressible transfer term appears in the transport equation of the com- 
pressible auto-correlation. As shown in Figure 9, the compressible transfer is positive for all spectral space. 
Hence, it is a term that is responsible for the production of compressible energy. 

The different contributors to T cc are plotted in Figure 10. It is clear that the two terms, T 2 CC and Tf c , 
are much larger than the others. Another term, 7\ cc , is a distant third in size. To illustrate the relative size of 
the smaller terms, we replotted these contributors in Figure 11 at the enlarged scale. Because two dominant 
terms have similar magnitudes but opposite signs, a strong cancelation between them is expected. Indeed 
for all cases considered, the term T$ c is always positive, whereas the other term T 5 CC is always negative. 
In fact, the cancelations are so complete, the summation of these two terms is now actually negligible in 
comparison with T x cc (Figure 12). The physical explanation for this 'almost perfect’ cancelation is that 
T 2 cc and Tf c are the terms that take into account the interactions between slowly varying incompressible 
modes and two compressible modes (namely the interactions between acoustic waves and a solenoidal field). 
This interaction results in the production of acoustic energy on the same wave-number but now in another 
propagation direction. For an isotropic redistribution of acoustic energy, this effect of reorientation does not 
affect the spectral distribution of energy and leads to a zero net balance. Consequently, the most important 
term in the transport equation of E cc is Tf c . Figure 13 demonstrates that this term is also much larger 
than all other contributors (T £ c , T£ c and Tjf c ). 
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Fig. 12. 

Comparing T§ S (K) to Tp c (K), it is clear that they have similar magnitude but with an opposite sign 
(Figure 14). These two terms are essentially responsible of the energy exchanges between the solenoidal and 
compressible parts. Specifically, Tp c is the “input” energy term on the compressible mode whereas T^ s is 
the “output” term in the equation of the solenoidal mode. Based on these results obtained in spectral space, 
we conclude that there is a local transfer of energy from the solenoidal mode to the compressible mode. This 
result will be further confirmed by our analysis in the second part of the paper. 

The total compressible transfer term T cc is dependent on the compressibility (see Figure 15). As 
expected, its magnitude increases with the increasing of turbulent Mach number. We note that there is 
a shift of the peak spectrum towards the large K. Based on the properties of T§ S {K), we expect that 
its compressible counterpart Tf :c (K) should also scale as . Since the spectra of divided by M t 
collapses (Figure 16), this term can not be responsible for the peak shift of T cc . Figure 17 showed that 
the ‘almost perfect’ cancelation between terms T? C {K) and Tg c (K) for all turbulent Mach numbers is 
considered. Figure 18 illustrates how the term T% C {K) depends on M t . Figure 19 shows the dependence 
of Tf c {K) and Tf c (K) on the values of turbulent Mach number. Although these two terms maintain 
opposite signs for all spectral space, the magnitudes and shape of T^ C (K) and T§ c {K) are clearly different. 
The ‘imperfect’ cancelation between these two terms leads to a cascade type of compressible energy transfer 
(Figure 20). These two terms involve ( E cc ) 2 and are important contributors at high turbulent Mach 
numbers. The interaction among the compressible mode begins to have influence. This cascade mechanism 
will be investigated in the next section. 
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Fig. 13. 


4. Triadic Interactions. The most fundamental building block of the energy transfer process is the 
triadic interactions. Specifically, we are interested in the energy transfer for a given mode K due to its 
interactions with all the pairs of modes P and Q = K - P that form a triangle with K. For this reason, we 
introduce the triadic energy transfer function, T(K , P, Q), according to 


(27) T SS {K)= 5Z T SS (K, P, Q),and 

P,Q= |K-P| 

(28) T CC (K) = J2 T cc (K,P,Q). 

P,Q=IK-P] 

Here T(K> P, Q) is defined as energy transfer to K due to triads with one leg in Q and the other in P. The 
average procedure is performed over spherical shell since the turbulence is isotropic. 

An examination of the purely incompressible contributors (T§ S {K, P, Q)) reproduces the results of in- 
compressible turbulence (Domaradzki and Rogallo (1990), Yeung et al. (1991, 1995), Ohkitani and Kida 
(1992), Zhou (1993a-b) and Zhou et al. (1996)) and again indicates that the purely solenoidal triadic energy 
transfer is not affected by compressible effects. The triadic solenoidal transfer T SS (K , P, Q) (with the com- 
pressible terms included) is essentially the same for a wide range of turbulent Mach number values (Figure 
21). Although Figure 21 is only for P = 512 and Q = 128, we have examined other values of P and Q and 
found that our conclusion does not change with turbulent Mach number. As a result, the compressibility 
has very little influence on the solenoidal triadic interactions. Figure 22 is a typical plot showing how the 
structure of T SS (K, P, Q) changes with various Q values (P = 512, M t = 10 -2 ). Again, this result is the 
same as that of incompressible turbulence. 


12 




— ■*— T ss 

1 5 

f - *p l * ' 


3 10 s 



We have found from the previous section that particular attention should be paid to the T§ s term since 
it is the term that is responsible for ‘output’ energy from solenoidal to compressible mode. As the Mach 
number increases, the magnitude of the ‘output’ energy increases but the basic structure remains (Figure 
23). Since this term represents an energy output at a localized spectral region, we refer to it as the radiative 
(emission) triadic energy transfer. Figure 24 illustrates this type of interaction from various Q values at 
the given value of P (P - 512) for M t = 10~ 2 . It is clear that the triadic interaction of this term is quite 
different from those of purely incompressible terms. 

We now turn our attention to the triadic interactions in compressible energy transfer, T CC (K,P,Q). 
In Figure 25, we present T CC (K,P,Q) for various Q values when P is in the inertial range (P = 512). 
For this low Mach number lO" 2 , we observe that the structures of T CC (K,P,Q ) are rather similar for 
differing Q values. All of them show the radiative (absorption) type of energy transfer. Figure 26 shows 
that Tf c (K, P, Q) is the dominant contributor to the compressible triadic energy transfer function. In 
fact., comparing Figures 24 and 26, we found that the absorption types of triadic energy transfer functions 
have the same magnitude but opposite sign as those of the emission type {T$ s term). Figure 27 further 
demonstrates that the triadic interactions Tf c {K,P,Q) ~ -Tg s (K,P,Q) for all turbulent Mach numbers 
are under consideration. As a result, we conclude that all compressible energy has been transferred locally 
(in spectral space) from the solenoidal component. 

From the previous section, we have found that the sum of Tf c (K, P, Q)+Tf c (K, P, Q ) is an important 
contributor to the compressible energy transfer. At low Mach number (M, = 10 -2 ), terms Tf c (K,P,Q) + 
Tg c {K, P, Q) are small and their triadic energy transfer terms show only a very weak energy cascade (Figure 


13 





Figure 1 5: Compressible transfer term at different N4 


Fig. 15. 


28). However, this situation changes rapidly as the Mach number increases. Indeed, the compressible energy 
cascade can be seen in Figures 29 and 30 where Tf c (K, P, Q) + Tg c (K, P, Q) are plotted for several higher 
Mach numbers. From this analysis, we conclude that at high Mach number the cascade of compressible tur- 
bulence is a direct result of the fact Tf c (K, P, Q) + Tg c {K, P, Q) > Tf c (K, P, Q). To further demonstrate 
this point, we plot the total compressible energy transfer term, T CC (K, P, Q), at different Mach numbers. It 
is clear that T CC (K , P, Q) changes its characteristic features from radiative to cascade as the turbulent Mach 
number increases (Figures 31 and 32). This is a result that can not be observed from the total compressible 
energy transfer function T CC {K). 
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Figure 22:Triadic soienoidal transfer term for different Q 
Fig. 23. 
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Fig. 31. 
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5. Conclusion. We have investigated the energy transfer process of compressible turbulence using a 
two-point closure theory, Eddy-Damped-Quasi-Normal-Markovian (EDQNM) closure, as well as the Helmholtz 
decomposition, a method that separates the compressible and solenoidal modes. We focused on the following 
issues: (1) What is the mechanism of energy exchange between the solenoidal and compressible modes, and 
(2) Is there an energy cascade in the compressible energy transfer process? We found that the compressible 
energy is transferred locally from the solenoidal part to the compressible part. We also found that there is 
an energy cascade of the compressible mode for high turbulent Mach number (M t > 0.5). Since we assume 
that the compressibility is weak, the magnitude of the compressible (radiative or cascade) transfer is much 
smaller than that of solenoidal cascade. We confirmed these results by studying the triadic energy transfer 
function, a most fundamental aspect of the energy transfer. 
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